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Empirically, Type la supernovae are the most useful, precise, and mature tools for determin- 
ing astronomical distances. Acting as calibrated candles they revealed the presence of dark 
energy and are being used to measure its properties. However, the nature of the Type la 
explosion, and the progenitors involved, have remained elusive, even after seven decades of 
research. But now new large surveys are bringing about a paradigm shift — we can finally 
compare samples of hundreds of supernovae to isolate critical variables. As a result of this, 
and advances in modeling, breakthroughs in understanding all aspects of these supernovae 
are finally starting to happen. 

"Guest stars" (who could have imagined they were distant stellar explosions) have been sur- 
prising humans for at least 950 years, but probably far longer. They amazed and confounded the 
likes of Tycho, Kepler, and Galileo, to name a few. But it was not until the separation of these 
events into novae and supernovae by Baade and Zwicky that progress understanding them began 
in earnest^. This process of splitting a diverse group into related subsamples to yield insights into 
their origin would be repeated again and again over the years, first by Minkowski when he sepa- 
rated supernovae of Type I (no hydrogen in their spectra) from Type II (have hydrogen) 2 , and then 
by Elias et al. when they determined that Type la supernovae (SNe la) were distinct 3 . We now 
define Type la supernovae as those without hydrogen or helium in their spectra, but with strong 
Sill, which has observed absorption lines at 6150, 5800 and 4000A 4 . The progenitor of a Type la 
supernova has never been seen, but the consensus is that they are the result of the thermonuclear 
explosion of a degenerate carbon-oxygen white dwarf (WD) star 5 gaining mass in a binary sys- 
tem^, in contrast to the other types of supernovae, which are thought to result from the collapse of 
the core of a massive star. Much about them remains unknown — the process of categorizing and 
trying to understand them continues. 

As distance indicators, SNe la are capable of better than 6% precision in distance after em- 
pirical corrections are made for their lightcurve shape and color. Just over a decade ago, they were 
used to deliver perhaps the most shocking scientific discovery in the last half-century: the expan- 
sion of the universe is not slowing, but is instead accelerating 7 8 . The "dark energy" required to 
power this acceleration is prodigious, and indeed must be the dominant constituent of the universe 
by a factor of 17 over normal matter, and 3 over dark matter (which, despite the name, is probably 
unrelated). But because dark energy is so dilute (one cubic meter contains the energy equivalent to 
a few atoms of hydrogen), its effects can only be seen on vast scales. We must compare supernovae 
across ^10 Gyr of cosmic time, during which stars and galaxies have undergone substantial evo- 
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lution. Therefore it is critical to understand supernovae if they are to be relied on as the primary 
distance probe in astrophysics. 

We bein by reviewing SN la lightcurves and spectra (^T]) and their use in cosmology 
This has been refined to the point that SNe la are now sensitive to subtle, systematic effects 
However, many of these are being addressed by comparing large subsamples of hundreds of Sis 
la, possible for the first time only in the last few years. 

New multidimensional simulations are revealing once-elusive details of the explosion, such 
as how asymmetry contributes to the diversity seen in SNe la (Q. Advances in understanding 
progenitors have been proceeding apace, led by new observations. The properties of SNe la are 
being correlated with the stellar populations in their host galaxies, yielding real progenitor con- 
straints (Q. Dramatic new discoveries, like supernovae that appear to require a progenitor above 
the Chandrasekhar mass, are challenging pre-existing theories. Evidence has been mounting that 
the long-favored scenario of accretion from a main sequence or red giant companion star cannot 
account for all SNe la. Instead, the white dwarf merger and even the sub-Chandrasekhar mass 
explosion models have been experiencing a revival. 

Finally, we conclude that empirical SN la cosmology remains solid, but is being continually 
refined by taking into account ever-more- subtle effects and new knowledge. But our understanding 
the the progenitors of SNe la and the explosion process may be undergoing an evolution, as the 
data from vast new surveys with thousands of supernovae are becoming available (@. 

1 Lightcurves and spectra 

In a SN la, a wave of thermonuclear fusion rips through a degenerate white dwarf star, synthesizing 
iron-peak elements (Ni, Co, Fe) in the dense inner regions, intermediate mass elements (Si, S, Ca, 
Mg, O) where burning is incomplete, and sometimes leaving unburned material (C,0) near the 
outer layers^. Though the explosion provides the kinetic energy of the SN la, and unbinds the 
white dwarf, this is not what we see as the supernova. The lightcurve is powered by the radioactive 
decay of 56 Ni (half-life 6.1 days) to 56 Co, and ultimately to 56 Fe (half-life 77 days). Gamma rays 
produced in the decays are thermalized, and at peak light ~ 85% of the light output of the SN la is 
in the optical, peaking at 4000A, with the remainder mostly radiated in the near-ultraviolet (NUV) 
and near- infrared (NIR) 1 . 

Figure [T] shows the quasi-bolometric (UV-optical-IR) lightcurve evolution of a typical SN la, 
SN 2003du^Talong with representative spectra. SNe la rise to maximum light over ~ 13 — 23 days 
in the £?-band, with a typical rise time of about 17.5 days 12 . Initially, the time for photons to diffuse 
out of the dense ejecta is high, so that the rate of deposition of energy by radioactive decay exceeds 
the energy radiated by the SN. Thus, at early times, spectra show only absorption lines probing the 
outer layers of the SN. As the SN expands and thins, the light-emitting region (the photosphere) 
recedes (in mass or velocity coordinates), and spectra start to probe deeper layers. Eventually the 
ejecta reach a point where the rate of energy deposition by radioactive decay is equal to the radiated 
luminosity of the SN, and the SN starts to decline in luminosity 13 . Around or after maximum light, 
spectra show P-Cygni profiles - emission at the rest wavelength (initially weak), and blueshifted 
absorption. Weeks after maximum light, the spectra start becoming dominated by scattering from 
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permitted lines (at least in the blue) 14 and ultimately become dominated by emission features as 
the ejecta start to become optically thin, and the SN makes the transition from the photospheric to 
the nebular phase. At late times (hundreds of days), gamma rays can freely escape, but positrons 
may be trapped^, and the lightcurve slope may match the slope expected from the decay of 56 Co, 
or can be steeper than this in the case of incomplete trapping^. Late time spectra are dominated by 
emission lines from iron-peak elements, synthesized in the deepest regions where the white dwarf 
was densest. 

2 Cosmology 

Dark energy is often characterized by its ratio of pressure to density, i.e., its equation of state, 
w — P/pc 2 (Fig. g). If a is an arbitrary length scale in the universe, then the density of some 
component of the universe evolves as p oc a~ 3 ^ w \ Normal matter has w = 0; it dilutes with 
volume as space expands. If dark energy has an equation of state w — — 1, corresponding to 
Einstein's cosmological constant, A, then it has the strange property that the energy density does 
not dilute as the universe expands; it must be a property of the vacuum. Other possibilities exist: 
if w < —1 the dark energy density is ever-growing, and will ultimately result in a "big rip," 
destroying galaxies and even subatomic particles. However, w can also be a scalar field with 
w > — 1 (generally known as quintessence), and the value of w can evolve over time. 

Measurements of (w) or w(z) are obtained by building a map of the history of the expansion 
of the universe using SNe la as standardized candles. SNe la can show a factor of 10 or more 
difference in peak luminosity, but the luminosity is correlated with the time it takes the supernova to 
rise and fall in brightness 17 . Therefore, the width of the lightcurve is measured and used to correct 
the peak luminosity. Common parametrizations include the "stretch factor 8 ", s, proportional to 
lightcurve width, or Ami 5 (5), the number of magnitudes the S-band lightcurve falls in 15 days 
after peak brightness, inversely proportional to lightcurve width 17 . A correction must also be made 
for color, since redder SNe are dimmer, both intrinsically, and due to dust ^3. Various techniques 
are used to determine these parameters by fitting lightcurve models to the data, but leading fitters 
include MLCS2k2^, SALT2 E H, and SiFTCP. 

To determine cosmological parameters, an observed Hubble diagram (distance versus red- 
shift) is constructed, and cosmological parameters are varied in a model that is fitted to the data. If 
we express distances as magnitudes, as a distance estimator, we can use^S; 

I^b = m* B - M + a(s - 1) - /3c, (1) 

where m* B is the peak magnitude of the SN in the £?-band (the blue filter where the SN is brightest), 
s is the lightcurve stretch, c is t he co lor of the SN (a linear combination of U — B and B — V, 
relative to some reference color 1 ^^), a is the slope of the stretch-luminosity relation, f3 is the 
slope of the color-luminosity relation, and M is a measure of the absolute magnitude of the SN 
combined with the Hubble constant. Note that since relative magnitudes are used, neither the 
absolute magnitude of the SN, nor the value of the Hubble constant must be known. m* B , s 9 and c 
are measured from a fit to each SN lightcurve, while M, a, and f3 are constants determined from 
the overall cosmological fit. 
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3 Systematic uncertainties 

Cosmological studies using SNe la, have reached a point where over most redshift ranges system- 
atic errors (those that affect many measurements simultaneously in a correlated way) dominate 
statistical errors (i.e. those that are reduced by y^]v)|24|25|26| — see pjg |^J Many systematic un- 
certainties can be lowered with new methods, improved statistics, by comparing subsamples of 
supernovae, or breakthroughs in understanding. The primary systematic uncertainties affecting 
SNe la are survey-dependent, but there is general consensus that the largest are: calibration to 
the historic Landolt photometric system, treatment of the ultraviolet, reddening due to dust, the 
differences supernovae show with respect to environment, and the possible evolution of SNe with 
redshift. See Table 1 for a summary of recent SN la constraints on (w) and dominant systematic 
uncertainties. 

Calibration As supernovae are observed at high redshift, the region of the spectrum seen through 
a given broadband filter changes. A correction must be made for this, the k-correction, which 
requires knowledge of the spectral energy distribution of an average supernova^, and precise 
knowledge of the filter transmission curves and calibration system used for each supernova. Un- 
fortunately, many historical low-redshift supernovae were transformed to the Landolt photometric 
standard system, which is poorly understood and no longer reproducible, as the filter and instru- 
ment transmission curves used to establish it are not known to modern precision and no longer 
exist. The fact that the low-z SN sample is on the Landolt system has forced many high redshift 
surveys to transform to it, incurring a systematic uncertainty in zero points and color terms^S. 
These systematic uncertainties will be reduced once a large new low-redshift sample of SNe is 
assembled, calibrated onto a new system (e.g. Sloan). 

The Ultraviolet At redshifts above z ~ 0.2, the restframe [/-band and ultraviolet is redshifted 
into the optical region of the spectrum. Unfortunately, there is a great deal of scatter in restframe 
U observations. Some surveys rely heavily on the restframe U, so uncertainties from this region 
can translate into dominant systematic uncertainties 24 . In fact, differences seen between the dis- 
tances determined by different lightcurve fitters, e.g MLCS2k2 l ^l, SALT2 E 3, and SiFTCP boil 
down largely to the way they treat the restframe U -band. 

There are two possibilities for [/-band uncertainty: (1) SNe intrinsically show more varia- 
tion in [Z 28 , and/or (2) the [/-band is poorly calibrated. The Supernova Legacy Survey (SNLS) 
has shown that when longer- wavelength high redshift observations from a well-calibrated single 
survey are k-corrected to restframe U -band, the dispersion is more than 3 times lower than for low- 
redshift data where the U -band is observed directly 23 . Therefore, discrepancies are most likely due 
to problematic low redshift observations where many factors conspire to make U observations dif- 
ficult, including atmospheric variation, extinction, nonstandard filters, and poor calibrators 29 . This 
systematic can be dealt with by building a new, better-calibrated low redshift sample, and possibly 
by avoiding the restframe U -band. 

Reddening Lightcurve fitters must correct for the fact that redder supernovae are dimmer. This 
is due to a combination of an intrinsic color-luminosity relation (faint supernovae are intrinsically 
redCED), anc [ reddening due to dust. While there is consensus that the two effects ought to be cor- 
rected for independently, there is disagreement over whether this can be practically achieved given 



4 



current data limitations. 

MLCS2k2 attempts to separate intrinsic and dust reddening and make different corrections 
for eactP2. SALT2 and SiFTO work under the assumption that it is difficult or impossible to 
separate the effects of dust and intrinsic reddening, so they make an empirical reddening correction 
by solving for the slope of the color- luminosity relation, f3. If this reddening were only due to dust, 
then /3 = Rb = Rv + 1- For Milky Way dust, R B = 4.1. However, as shown in Fig. [3} shallower 
slopes are found for /?. This may result from the conflation of dust and intrinsic SN reddening, and 
would indicate that the intrinsic SN la reddening-luminosity relation has a shallower slope than 
the dust relation. However, since the lightcurve shape-luminosity relation is already factored out 
in this method^, it means there must be a component of intrinsic color that does not correlate with 
lightcurve shape 19 . An alternative is that the dust along the line of sight to SNe la is intrinsically 
different, or that scattering effects in the circumstellar environment result in a different apparent 
reddening law^. 

One way to separate the intrinsic/dust reddening degeneracy is through infrared observations, 
where the effects of dust are minimized. Infrared observations have shown that many SNe seem 
to have a sub-Milky Way value of R v , sometimes as low as R v = 1.5^21. There are hints that 
lightly reddened SNe may have an R v close to that in the Milky Way, while heavily extincted SNe 
la have lower values (Fig. 3)^. However, this is opposite to earlier findings from optical data^, 
and the small sample size makes it difficult to draw reliable conclusions. 

There are indications that supernovae with higher velocity ejecta have lower inferred R v 
values 35 (Fig. [3]). Since the color is correlated with physical SN features, this may indicate that 
sometimes an intrinsic color difference in SNe is wrongly ascribed to dust reddening (though this 
can be mitigated by making a color cut), or perhaps that there are different progenitor scenarios, 
with different circumstellar dust scattering properties, that produce SNe with different velocities. 

Various authors have exploited the low sensitivity to dust in the IR to make infrared Hubble 
diagrams. This is doubly attractive, because SNe la intrinsically show reduced scatter in the IR^ 
— in the H band SNe la have an rms dispersion of only 0.15 mag without correction for lightcurve 
shape or color 37 . This is theoretically expected, in part because less luminous SNe are cooler and 
radiate a larger fraction of their luminosity in the IR^. While heroic attempts have been made 
to produce high redshift IR Hubble diagrams of SNe la 39 , the massive amount of observing time 
required, and the fact that the restframe IR is redshifted to even longer wavelengths, has made 
serious cosmological constraints from this method elusive. 

Populations and evolution Even before SNe la were used for cosmology, it was known that the 
most luminous SNe la (those with the broadest lightcurves), occur only in late- type galaxies^. 
Likewise, subluminous SNe la are preferentially found in galaxies with a significant old popula- 
tion^. SN lightcurve width and luminosity have now been shown to correlate with host galaxy star 
formation rate, galaxy mass, and metallicity, | 42|43|1Q | 44 l 

Since star formation increases by a factor of 10 up to redshift 1 .5, and bright, broad- lightcurve 
supernovae favor star forming hosts, it is expected that the mix of supernovae will change with red- 
shift. Indeed it has been shown that supernovae at z = 1 are intrinsically ~ 12% more luminous 
than local SNe la 45 . The spectra of high-redshift supernovae also show fewer intermediate mass 
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elements, consistent with the idea that they make more iron-peak elements to power their luminos- 
ity^. A changing mix of supernovae with redshift is not necessarily problematic for cosmology if 
lightcurve shape and color corrections allo w all supernovae to be corrected to the same absolute 
magnitude. Unfortunately, they do not^EH — as is shown in Fig. [4| supernovae in high and low 
mass galaxies each correct to an absolute magnitude different by 0.08 ± 0.02 mag. This trend is 
also present (though weaker), if SNe are split by host star formation rate, or metallicity. If not 
corrected, this can lead to a ^ 0.04 systematic error in ufiSl. Therefore, aside from color and 
lightcurve shape, a third correction, one for host galaxy properties, must be applied to SNe la to 
avoid systematic residuals with respect to the Hubble diagram. 

The physical or igins fo r the differences in supernova properties in spiral and elliptical galax- 
ies have been elusive^^H _ are th e y related to progenitor age, metallicity, or entirely different 
progenitors? The recently d iscov ered trend that galaxy mass may play the most significant role 
in determining SN properties'^^ would seem to implicate metallicity, since higher mass galaxies 
retain more metals in their deeper potential wells. Still, the age-metallicity degeneracy precludes 
firm conclusions on this evidence alone^. Theoretically, there are reasons that high metallic- 
ity progenitors should produce subluminous SNe Ia^ increased 22 Ne in high metallicity white 
dwarfs provides more neutrons during the nucleosynthesis that occurs during the explosion, and 
thus yields more stable 58 Ni and less of the radioactive lightcurve-powering 56 Ni. However, this 
has only a ~ 10% effect on SN luminosity, perhaps indicating that an as-yet unknown metallicity 
effect or progenitor age is the more important in controlling 56 Ni yield and SN luminosity^. 

If progenitor age is the predominant effect controlling SN luminosity, this has implications 
for SN progenitors. In the single degenerate (SD) scenario the time to explosion is controlled 
by the main sequence lifetime of the secondary (mass transferring) star, since the mass transfer 
timescale is thought to be negligible, only ~ 10 7 years. Therefore, in this model, the delay time 
is a direct indication of the secondary star's mass. In the double degenerate (DD) model, where 
two white dwarfs merge, gravitational wave radiation ultimately brings the stars together, but they 
may have a head start if they have drawn closer because at some time they have been orbiting in 
the common envelope of one of the stars as it is evolving. Age may still play a role in determining 
luminosity, because at early times the only white dwarfs will be massive ones derived from the 
more massive stars. 

4 Explosion theory and observational constraints 

In what might be called the standard model for a SN la, a CO white dwarf accretes matter until 
it compresses to the point that carbon is ignited just before the Chandrasekhar limit (note that the 
common misconception that the WD goes over the Chandrasekhar limit is wrong - this would 
lead to collapse to a neutron star). The evidence that the exploding star is a white dwarf is strong, 
albeit circumstantial: neither hydrogen nor helium is seen in the spectrum of a SN Ia^, SNe la 
can happen long after star formation has ceased, the explosive process may implicate degenerate 
matter, the energy obtained by the thermonuclear burning of a white dwarf minus the binding 
energy roughly matches the kinetic energy of SNe la, and s imulations of the process have been 
successful at reproducing SN la lightcurves and spectra'^EZE]. 

The white dwarf may have a "simmering" phase of order a thousand years following unstable 
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carbon ignition, where thermonuclear runaway is prevented by convection 53 . Ultimately, however, 
explosive burning is ignited and the white dwarf is incinerated in seconds. 

If a white dwarf near the Chandrasekhar mass is detonated (i.e. the burning occurs superson- 
ically), then the white dwarf burns at such a high density that the fusion products consist almost 
entirely of iron-peak elements^. This does not match the spectra or lightcurves of SNe la. A 
deflagration (subsonic burning), on the other hand, gives the SN time to pre-expand. Burning at 
a lower density can produce intermediate mass elements and reproduce many of the observational 
features of SNe la 51 . However, pure deflagrations fail to produce the high velocity material (Fig. [I]) 
that seems to be nearly ubiquitous in SNe la in the outer layers^, and in two- and three dimen- 
sional simulations they leave too much carbon and oxygen unburned^. Therefore, the consensus 
is that in a Chandrasekhar mass explosion, the flame must start out subsonically, but at some point 
become supersonic^. Though the physics remains poorly understood, models that start as a de- 
flagration but impose a transition to a detonation under certain conditions have been successful 
at reproducing normal SN la lightcurves and spectra, and even lightcurve width-luminosity rela- 
tions and metallicity effects 57 . Explosions dominated by deflagration produce more intermediate 
mass elements (dimmer SNe la), while those dominated by detonation produce brighter and more 
56 Ni-rich SNe la (Fig.]^ 3 . 

It is possible to gain insight into the explosion physics with spectropolarimetry or other tech- 
niques that reveal asphericity — recent work indicates that some of the dispersion in SN properties 
results from broken symmetries 57 . SN asymmetry can be measured via spectropolarimetry, since 
asymmetric electron scattering leads to p olariz ation vectors that do not cancel. Most normal SNe 
la are found to be spherically symmetric^^l The first convincing evidence for significant de- 
viations from spherical symmetry was seen in a subluminous SN Ia^, which may make sense if 
they are deflagration-dominated. However, supernovae with high velocity features often show even 
stronger spectropolarimetric signatures of asymmetry, possibly due to clumpy ejecta^^H 

Strangely, SNe la whose spectroscopic features start off with high velocity and evolve rapidly, 
often show nebular lines that appear redshifted, while SNe la with slower velocity evolution show 
blueshifted nebular lines 62 . This probably indicates asymmetry in the explosion, in qualitative 
agreement with models where a deflagration burns off-center, and is followed by a detonation. 

5 The progenitor question 

Even if there is agreement that the primary star in a SN la is a white dwarf, the identity of the 
secondary star in the binary system is uncertain. There are three broad classes of models: (1) 
single degenerate 6 , where the companion is a main sequence or red giant star that loses mass via 
either Roche lobe overflow, or a wind (symbiotic star), (2) double degenerate (DD), where two 
white dwarfs merge and explode'^SSD anc [ (3) sub-Chandra 65 , where a layer of helium builds up on 
the surface of a white dwarf below the Chandrasekhar mass until it detonates. 

Slow accretion In the single degenerate model, a white dwarf ignites carbon burning, causing it 
to explode, by accreting matter from a nondegenerate companion to near the Chandra mass. Since 
CO white dwarfs are thought to start no larger than ~ 1 — 1.2 solar masses, a significant amount of 
accretion from a secondary star is required. However hydrogen must be accreted at a rate between 
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10 -7 — 10~ 8 M yr" 1 to steadily burn to carbon or oxygen on the surface of the white dwari^. 
If it accretes faster, it will form a red giant-like envelope. Slower accretion is thought to lead to 
the build up of matter that results in novae, and mass being lost from the system. The apparently 
required steady burning generates supersoft X-rays, and as a result systems in this phase of their 
evolution are observable as supersoft X-ray sources 67 . However, the observed rates of sup ersoft 
sources in the Milky Way and external galaxies can only account for ~ 5% of the SN la rate 1 ^^, 
possibly favoring the double degenerate model, although many of these systems ought to pass 
through a supersoft phase as welp3 

White dwarf mergers In the double degenerate scenario, two white dwarfs merge to achieve 
the Chandra mass^MJ. The less massive white dwarf is disrupted into a disk that will eventually 
accrete onto the more massive white dwarf 7 ^. In this case, carbon may be ignited on the surface 
of the more massive white dwarf during the accretion process, resulting in non-explosive carbon 
burning, converting the star to an O-Ne white dwarf. This should ultimately undergo accretion 
induced collapse to a neutron star rather than producing a SN Ia^l, though it may be possible to 
avoid this fate under certain conditions 73 . The simulated merger of two equal mass 0.9 M WDs 
has been shown to lead to a SN la, although it was subluminous 74 , due to the low central density of 
a 0.9 M WD. More massive equal-mass mergers may produce normal, or even overluminous SNe 
la, but would be more rare. Lower-mass mergers should be common, though it is not clear if they 
lead to SNe la 75 . The merger of two white dwarfs is a very complicated three dimensional process 
(necessitating approximations in current models), and very few simulations have been completed. 

The WD merger scenario has a natural explanation for greater SN luminosity in young envi- 
ronments: younger, more massive stars produce more massive white dwarfs. Massive white dwarf 
mergers have more potential fuel than less massive mergers. 

Sub-Chandrasekhar mass explosions While approaching the Chandra mass is a convenient method 
for triggering carbon burning, there is no hard evidence that the Chandra mass is required. In the 
classical sub-Chandra models^, known as "double detonation" or "edge-lit detonations," a layer 
of accreted helium (~ O.2M ) is built up either by burning accreted hydrogen to helium or by 
accretion from a helium rich donor. When the pressure at the base of the helium layer reaches 
a critical threshold, it detonates, driving a shock into the core of the WD. This causes a second 
detonation, resulting in a flame propagating outward from the core (or near it), destroying the WD. 
Because the sub-Chandra white dwarf has a lower density throughout, a simple detonation does 
not burn the entire star to iron peak elements. This model thus has the advantage that an ad hoc 
deflagration to detonation transition is avoided. However, the significant outer helium layer is ef- 
ficiently burned to 56 Ni, resulting in early time spectra that should be rich in 56 Ni, which does not 
match observations^. If helium can detonate with a smaller layer, which some studies hint 77 , then 
sub-Chandra simulations can reproduce many of the observed properties of SNe la 78 . Some of the 
elegance of the model is undercut by the required narrow range of WD masses: ~ 0.9 — 1.1 M , 
in order to achieve the central densities necessary to produce iron-peak elements. 

Constraints on progenitors from rates Different progenitor scenarios lead to different delay 
times between the birth of the binary system and the explosion as a SN la 79 . So there is hope of 
constraining the progenitors by measuring the delay time distribution (DTD; Fig. [6]). One approach 
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is to measure the lag between the cosmic star formation rate and the SN la rate as a function of 
redshift. However, large uncertainties in each make this particularly difficult, leading authors to 
different conclusions even when using largely the same dataf^2IS3J 

Another approach is to use the relative SN rates in different types of galaxies to constrain the 
DTD. The fact that the SN la rate per unit mass is higher by a factor of ~ 20 in late-type galax- 
ies indicates that there is a significant population of short-lived thermonuclear SN progenitors 82 . 
While there has been much argument about whether the SN la DTD is bimodal^l, a leap forward 
came with the application of SED fitting techniques to SN host galaxies by modeling their stellar 
populations^. This has allowed more complicated DTDs to be determined: authors using likely 
SNe la from the SXDS^ (squares in Fig. [6]), or those in clusters 85 (filled circles), find a power 
law DTD proportional to t~ l . It is difficult to explain the power law DTD results using the single 
degenerate scenario alon e |79 | 84 i though there is at least one claim^. 

Host galaxy spectroscopy contains even more infor mation to constrain progenitors. New 
studies using different methods and different data sets^ES but the same galaxy fitting code, have 
determined the relative rates of supernovae in 0-0.42 Gy (prompt), 0.42 to 2.4 Gyr (medium), and 
>2.4 Gyr (delayed) bins (stars Fig [6]). They find that "prompt" and "delayed" supernovae are re- 
quired at several sigma. They also confirm earlier findings that bright, broad- lightcurve supernovae 
favor a prompt population, while dim, narrow-lightcurve SNe la favor a delayed population. 

Pinpointing the locations of SNe la relative to the stellar populations in which they reside 
can give some information on progenitors, even if it is imprecise due to the probability that the 
SN has migrated from its birth site. Comparing the distributions of SNe la to the blue light in 
galaxies reveals that even "prompt" SNe la are significantly delayed, with delay times 200-500 
MyjMl Additional constraints have been obtained by comparing SN la remnants in the Magellanic 
Clouds (triangles in Fig. [6]) to their local resolved stellar populations 90 . At a > 99% confidence 
level, "prompt" SNe la that explode within 330 Myr of star formation are required. 

Surviving or preceding material If white dwarfs grow to become SNe la via hydrogen accretion, 
then it ought to be possible to to see this hydrogen in the SN la in some form. However, despite 
decades of searching, no such smoking gun has been found. No SN la has ever been seen in the 
radio 91 , establishing an upper limit on steady mass loss of material before the SN of ~ 3 x 10 -8 
M yr _1 . This disfavors the symbiotic star hypothesis, where the white dwarf accretes matter from 
the stellar wind of its companion. If, on the other hand, white dwarf is accreting from Roche lobe 
overflow of a companion, the outer hydrogen layers of the secondary ought to become stripped 
and entrained in the SN ejecta, where it will show up at low velocities when the ejecta become 
optically thin 92 . However, in a couple of well observed SNe la, upper limits on the amount of 
hydrogen detected are 0.01 M , which would seem to rule out this scenario 50 , at least for these 
SNe. In a few cases, apparent SNe la have been seen interacting with pre-existing hydrogen, as in 
SN 2002ic^l, although it is not clear these are really SNe IaPI, and even if they are, these systems 
are the exception rather than the rule. 

The impact of the SN ejecta on any companion star would shock the ejecta, leaving observa- 
tional signatures along certain lines of sight for several days, long after the ejecta has overtaken the 
companion^. This would lead to distortions in the early-time lightcurves of SNe la around 10% 
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of the time, especially in the UV to B regions of the spectrum. If the companion star is a red giant, 
particularly large distortions ought to be seen - the fact that they never have been would seem to 
be problematic for this scenario^. 

Variable sodium lines have been seen in high resolution spectra of some SNe la, and their 
variability, assumed to be from ionization of circumstellar clouds by the SN, has been taken to 
be evidence of significant circumstellar material. This ostensibly favors the single degenerate 
scenario^, although maybe only in rare cases — variable sodium was only seen in 6% of cases in 
a sample of 31 observations, and then only in the most reddened SNe la 98 . Not all variable sodium 
SNe la are reddened, however, nor do all reddened supernovae show variable sodium^. There 
have also been doubts that this material is truly circumstellar^. 

SN remnants provide more clues to the origin of SNe la. They may reveal the progenitor's 
metallicity^, for example. But even more exciting is the prospect that they may house a surviving 
progenitor star. Tycho's remnant (SN 1572), long suspected of being a SN la, has been definitively 
revealed to be so by modern spectroscopic observations of the SN light, delayed in transit by 
reflections off of distant interstellar clouds^. In the single degenerate scenario, the secondary star 
should survive the SN explosion, but may be stripped of its outer layers, and heated, such that it 
takes of order 1000 years to return to thermal equilibrium^. One star, Tycho G, has been claimed 
as evidence for a surviving progenitor^, though there are many doubts about this hypothesis. It is 
too displaced from the center of the remnant, and its proper motion shows it did not come from that 
direction. Also, it was not apparently out of thermal equilibrium, and does not have the expected 
rotational characteristics^. 

Super-Chandra Beginning with SNLS-OSDSbbC 53 , aka SN 2003fg, several supernovae have shown 
luminosities indicating a Nickel yield alone of 1.3 - 1.8 M Q | 1Q6 | 1Q7 1 , possibly indicating a super- 
Chandrasekhar mass progenitor. Such energetic explosions ought to have more kinetic energy, and 
thus higher velocity ejecta, but this class of SNe counterintuitively have much slower than average 
ejecta. This may result from the increased binding energy of a super-Chandra WfP^, or from the 
ejecta plowing into a dense shell of circumstellar material resulting^. These SNe also show rare 
lines of CII near maximum light, possibly indicating less complete burning. They may also favor 
a young stellar population. This class of SNe la would seem to implicate the merger of two white 
dwarfs, but no simulation has yet reproduced one. 

6 Final Thoughts 

SN la are mature probes of cosmology: systematic uncertainties now dominate statistical uncer- 
tainties up to z = 1, and the systematic uncertainties are well characterized. Still, there is plenty 
of room for improvement in the use of SNe la as standard candles. 

It is possible that no single model can describe all SNe la. Perhaps some SNe la are sub- 
Chandra, some are the result of mergers, and some are caused by Roche Lobe overflow. The 
explosion of a white dwarf may be similar no matter how it is triggered. If it is not required that a 
model explain all SNe, this opens the door to scenarios that do not have high expected event rates. 

In the study of SNe la, we have left the "serendipity era," when we could only discover what 
nature happened to tell us through chance encounters with nearby supernovae, and we are entering 
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the "database era," where we can proactively ask and answer questions about supernovae by com- 
paring large subsamples. The legacy yields of the multi-year second-generation SN surveys are just 
becoming apparent. Meanwhile, some third-generation surveys like the Palomar Transient Factory, 
Pan-STARRS, the Dark Energy Survey, and Skymapper, have already started to produce results, 
and promise thousands more SNe la in the next few years. Finally, fourth-generation surveys, 
the Large Synoptic Survey Telescope (LSST), and possible space missions promise unimaginable 
riches on the horizon. New techniques like comparing large subsamples, using UV and IR data, 
state of the art calibration, using host galaxy observations, and 3d modeling, have started to yield 
results, but are still in their infancy. The next decade holds real promise of making serious progress 
understanding nearly every aspect of SNe la, from their explosion physics, to their progenitors, to 
their use as standard candles. And with this knowledge may come the key to unlocking the darkest 
secrets of dark energy. 
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Table 1 : Best-fit values of (w) and error estimates. For the CfA3/ESSENCE column, 
w is from Hicken et al. 2009^1, though uncertainties are from ESSENCE^, as they are 
stated to be similar but are not separately tabulated. The SDSS numbers 24 are for their 
SALT2 fit. Errors for each survey use their largest sample. 
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Phase relative to maximum (days) Restframe Wavelength [A] 

Figure 1: UVOIR lightcurve and selected spectra from the normal Type la supernova 
SN 2003du. (a) Quasi-bolometric (UV-optical-infrared) lightcurve for SN 2003du E3] . Or- 
ange points mark the dates of spectra shown in panel b. (b) Spectra of SN 2003du [H] , 
with phases relative to maximum light listed in orange. Characteristic intermediate mass 
element features seen in maximum light spectra of SNe la are noted in green: Call, SN, 
and Sill. Before maximum light, when we are seeing the outermost layers of the ejecta, 
Call is seen at high velocity, a common feature in SNe la. SN ejecta thin with time, and 
we see to deeper layers. As SNe leave the photospheric phase and enter the nebular 
phase, emission features start to dominate. 
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Figure 2: Improvements in dark energy equation of state constraints from Type la 
supernovae over time. The top two panels show the remarkable improvements in w 
made using SNe la over the past decade (assuming a flat universe), the bottom two show 
the importance of improving systematics in the early years of the next decade 111 , (a) 
fl M - w statistical-only constraints circa 1998^. (b) By 2008, Kowalski et al., combining 
many data sets, showed that systematic uncertainties are significant 113 , (c) Expected 
constraints for the year 5 results of SNLS, assuming additional low-z SNe, and double 
the number of z > 1 SNe from HST, assuming there is no improvement in systematic 
uncertainties from the 3rd year result, (d) Assumes the low-z data are on the Sloan 
photometric system, and a factor of two improvement in measurements of fundamental 
flux standards. Judged by the area of the inner 68.3% contour, the improvement from the 
1998 results is a factor of 3, 5, and 10, including systematics. 



14 



I — ■ — I 



F10avg. w/o 06X & 05A 




NG08 E(B-V) < 0.7 
NG08 E(B-V) < 0.25 
W09 Normal Velocity 
W09 High Velocity 
L10 Star Forming 
UOLowSF 
S10 Star Forming 
S10 LowSF 
• 1 



F10SN 2006X 
F10SN 2005A 
E06 SN 2003cg 
E08 SN 2002cv 
K07SN 2001 el 
K06SN 1999cl 



1.0 1.5 2.0 2.5 

R v or |3 - 1 



3.0 



3.5 



4.0 



Figure 3: Determination of the reddening law along the line of site to SNe la. In most 
cases R v = A V /E(B - V) is shown, but in the S1 data, is shown, equivalent to R v if 
the colors of supernovae after lightcurve shape correction are due to dust reddening. Dif- 
ferent symbols/colors correspond to different techniques. The most reliable is using u to 
IR data to determine the reddening law (red stars), but this can only be done for a few in- 
dividual SNe la. Data is from K06E3, K07E21, E08E2, and F10 33 . The average from 
the Carnegie Supernova Program (F10) IR data, excluding their most reddened events, 
is shown as a red circle, the only determination consistent with Milky Way dust (dotted 
line). However, NG08 34 (green triangles) come to the opposite conclusion (albeit using 
only UBVRI data): when they restrict themselves to the least reddened supernovae, they 
find a lower value of R v . Blue squares show reddening behavior determined by solving 
for the reddening law for ensembles of SNe that would minimize their residual on the Hub- 
ble diagram. Results split by star forming galaxies and those with low star formation are 
shown from the SDSS-II SN Survey (LIO 011 ) and the SNLS (S10 44 ). There are some indi- 
cations that solving for reddening by minimizing Hubble residuals is incompatible with IR 
determinations 33 . Finally, W09ESfind that supernovae with high velocity features in their 
spectra have a lower R v , perhaps an indication that there is a component of reddening 
that is intrinsic to the SN. 
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Figure 4: Residuals from the Hubble diagram for Supernova Legacy Survey Type la 
supernovae as a function of host galaxy mass after correction for lightcurve shape 
and color. Figure from Sullivan et al. 44 Red points are averages. Histograms on the right 
show the differences for high and low mass host galaxies, split at M = 10 10 M . Horizontal 
lines show the average residual for each histogram. SNe la in high and low mass hosts 
correct to an absolute magnitude different by 0.08 ± 0.02 mag (4a). Supernovae in low 
mass galaxies are on average brighter than those in high mass galaxies before correction, 
but are dimmer after s and c correction. This galaxy-dependent residual can be corrected 
by taking host galaxy information into account. 
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Figure 5: Two dimensional Type la supernova models showing the relative impor- 
tance of deflagrations, detonations, and asymmetry. Figure from Kasen et al. 57 Blue 
shows intermediate mass elements (Si, S, Ca), green is stable iron group elements, and 
red is 56 Ni. An initial deflagration wave produces turbulent instabilities, but a later deto- 
nation wave burns much of the remaining fuel to 56 Ni. In low density outer regions the 
detonation produces only intermediate mass elements. Weak detonations produce the 
lower 56 Ni mass, less luminous SNe la. 
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Figure 6: Delay time distribution for Type la supernovae. Adapted from Maoz, Sharon, 
and Gal-Yam^. All solid circles except the first are delay times derived from cluster SN 
rates, assuming formation at z = 3. The first solid circle results from a constraint on the 
observed iron-to-stellar mass ratio in clusters. Blue stars are derived from SN rates and 
galaxy stellar populations in the Lick Observatory Supernova Survey^ 1 . Green triangles 
(including an arrow denoting the 95% confidence upper limit) are constraints from SN 
remnants in the Magellanic Clouds. Red squares are based on SN la candidates in E/SO 
galaxies at z = 0.4 - 1.2 84 . Curves are power laws t~ 1A (solid), and t -1 - 3 (dashed), 
constrained to pass through the last point. The prediction 79 from DD mergers is t~ l 
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Figure 7: Schematic illustration of the arguments for and against a given progenitor 
scenario being the predominant mechanism for SNe la. Arguments on the left strongly 
disfavor a certain scenario, while arguments on the right strongly favor that model. Note 
that an argument against a given scenario need not be an equally strong argument for 
another scenario (i.e. the diagram need not be symmetric). Though all evidence is subject 
to theoretical interpretation, italics indicate exclusively theoretical arguments. Relative 
rankings reflect an effort to distill community consensus, but are subjective. Arguments 
and references are described in the text. "Population effects" refers to the fact that more 
luminous SNe la are found in younger environments. "High M WDs seen" refers to the 
fact that CO WDs > 1.2 M exist in nature, thought to only be possible by slow accretion. 
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